Quasi-periodic spin chains in a magnetic field 
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We study the interplay between a (quasi) periodic coupling array and an external magnetic field 
in a spin-| XXZ chain. A new class of magnetization plateaux are obtained by means of Abelian 
bosonization methods which give rise to a sufficient quantization condition. The investigation of 
magnetic phase diagrams via exact diagonalization of finite clusters finds a complete agreement with 
the continuum treatment in a variety of situations. 

PACS numbers: 7f.10.Fd, 71.10.Pm, 75.60.Ej. 
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The magnetic properties of quasi-crystals have become 
a fundamental issue of study since their discovery in 1984 
[El. A variety of theoretical efforts, ranging from renor- 
malization group (RG) analysis of Ising models in Pen- 
rose lattices to exact solutions of both Ising and XY 
Fibonacci spin chains HH have revealed fairly intricate 
magnetic orderings associated to the quasi-periodicity of 
these structures. The non-metallic spin exchange mech- 
anism implicit in those studies has been evidenced in re- 
cently synthesized rare earth (TV) ZnMg-7?. quasi-crystals 
(see e.g. Q) whose TZ elements have well localized 4/ 
magnetic moments. 

Bolstered by these latter findings and as a further step 
within the line of the local moment descriptions referred 
to above, here we consider the ordering of quasi-periodic 
spin-i XXZ chains in a magnetic field to elucidate the 
quantization conditions of massive spin excitations or 
magnetization plateaux. In periodic systems, this issue 
has received systematic attention in the last few years 
from both experimental and theoretical points of view 
(see e.g. 0). In this letter, we are specifically interested 
in studying the antiferromagnetic system 
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denote the spin-^ matrices involved in 



where S X ,S V ,S 
the standard XXZ Hamiltonian (e„ — ) in a magnetic 
field h applied along the anisotropy direction ( |A| < 1 ). 
Here, the coupling modulation is introduced via the 
e„ parameters defined as e n = Ylv ^ cos (2 71 " ) , so 
quasi-periodicity arises upon choosing an irrational sub- 
set of frequencies w„ with amplitudes 8 V . 

The interest of (Q) stems partly from the widespread 
applications of Id Hamiltonians in the description of ar- 
tificially grown quasi-periodic heterostructuresjKj, quan- 
tum dot crystals |J and magnetic multilayers ]1Q] . Also, 
recent investigations of quasi-periodicity involving cither 
the couplings jll| or the magnetic field jl2| , have been ad- 
dressed using Abelian bosonization along with RG and 



numerical techniques. Here we focus on the combined 
effect of a quasi-periodic exchange modulation under a 
uniform magnetic field. 

Of particular importance are the rational frequencies 
of ([!]), not only as a way to approach the quasi-periodic 
limit, but also because they allow for a thorough numeri- 
cal verification of a novel situation (see Jl3|-p^[ for related 
work). As we shall see, although the allowed fractional 
plateaux predicted in the present case fall into the classi- 
fication provided by the generalized Lieb-Schultz-Mattis 
theorem [^6|, a bosonization approach to (|l|) yields an 
alternative scenario not envisaged in previous studies 
Jl^,^8|. This will be reflected in the appearance of mag- 
netization plateaux associated to each of the frequencies 
present in ([j]) . To strengthen the potential interest of our 
results, we show how a simple two frequency model ex- 
hibits a magnetization curve with two wide plateaux at 
1/4 and 3/4 of saturation, a situation which is highly 
reminiscent of that observed in magnetization experi- 
ments on NH 4 CuCl 3 0. 

Following the standard bosonization procedure (see 
e.g. p0[|), the continuum limit of the XXZ Hamiltonian 
in the presence of an external magnetic field h is given 
by the Tomonaga-Luttinger Hamiltonian 
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The bosonic field cj> and its dual 4> are given by the sum 
and difference of the light-cone components, respectively. 
The constant K = K((M),A) governs the conformal 
dimensions of the bosonic vertex operators and can be 
obtained exactly from the Bethe Ansatz solution of the 
XXZ chain (sec e.g. for a detailed summary). One 
has K = 1 for the SU(2) symmetric case (A = 1) and 
it is related to the radius R of Q by R- 1 = 2ttR 2 . In 
terms of these fields, the spin operators read 
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where the colons denote normal ordering with respect to 
the ground state with magnetization (M). The Fermi 
momentum kp is related to the magnetization of the 
chain as kp = (1— (M))ir/2. Either the XXZ anisotropy 
or an external magnetic field modify the scaling dimen- 
sions of the physical fields through K and the commen- 
surability properties of the spin operators, as can be seen 
from (^), (|). The non- universal constants a, b and c can 
be in general computed numerically (see e.g. for the 
case of zero magnetic field) and in particular the constant 
b has been obtained exactly in [^2| . 

For the sake of clarity, let us first consider the effect 
introduced by a single frequency term in Eq. (fil), i.e. e n — 
6 cos (2tt uj n) . Thus, using Eqs. (|) and (|) it follows 
that the relevant part of the interaction term in H qp = 
H + H int reads 

Him = ^2 C ° S ( 2lTU;X ) X l(9x4>) 2 + ^2{d x 4>) 2 

X 

+ A 3 cos(2fc F a; + %/2tt0) + A 4 sin(2fc F a: + V2tt</>) (5) 

where Aj oc # , i = 1, •••,4. 

As in previous analysis [p~6|— p~8| one can readily obtain 
the necessary quantization condition for the appearance 
of a plateau by looking at the commensurability of the 
relevant operators. In the present case we need only to 
consider the vertex operators exp ±i\Z2n(j) of scaling di- 
mension d = l/(47ri? 2 ). Therefore, we obtain that (M) 
should satisfy 



(M) = ±(2w - 1) , 



(G) 



in order for a plateau to be present. The fulfillment of 
this condition opens a spin gap excitation since the op- 
erator in question is relevant at least for < A < 1. In 
the region — 1 < A < a critical curve appears which 
can be determined from the Bethe-Ansatz solution for 
R((M), A) [fL8| . Furthermore, the gap width can be eas- 
ily computed to scale as ^A 2- "*). Notice that the per- 
turbations Ai,2 do not play a role here since they are 
incommensurate whenever (Q) holds. 

In fact, these expectations turn out to apply well above 
the weak coupling regime discussed so far. In Fig. 1 we 
show the magnetization phase diagram resulting from ex- 
act diagonalization of fairly large XX chains (10 5 spins) 
by setting uj = 13/21 through a wide range of couplings. 
Notice that for rational frequencies, ui G Q and S — > 5 C , 
for some appropriate value of S c which depends on uj, the 
chain breaks up into a periodic collection of finite seg- 
ments which naturally yield additional plateaux of ra- 
tional values. Being the unit cell composed of 21 spins 
one would naively expect plateaux to appear for values of 
M = (2n+ 1)/21, n = 0, • • • , 10, but as can be seen from 
Fig. 1 this is not the case due to the non-trivial structure 
of the unit cell and only some of these values are indeed 
observed. Interestingly, Eq. (^) yet remains robust all the 
way through the decoupling point and the corresponding 
plateau corresponds to the most prominent. 



These observations were also corroborated on XXZ 
chains. However, owing to the large spaces involved in 
their diagonalization now we content ourselves with mod- 
erate lengths L (up to 24 spins). After resorting to the 
Lanczos method [ p3[ in each of the magnetization sub- 
spaces with S z € {0, 1, • • • , L/2} , we built up the mag- 
netization contours shown in Fig. 2 using uj = 5/8 for 
L = 8, 16, 24 . As expected, the massless XXZ excita- 
tions around (M) = 1/4 render size effects rather notice- 
able within small coupling regions. Nevertheless, already 
for 5 > 0.2 they clearly become less pronounced, thus 
lending our results further support to the bosonization 
picture. 

A word of caution should be added here, namely, the 
importance of periodic boundary conditions (PBC) in 
testing the analytic approach via exact diagonalization 
of small systems. Already at the level of a simple dimer- 
ized chain (uj = 1/2), PBC become crucial. In fact, 
the numerical analysis of this latter situation using open 
boundary conditions shows that the well known (M ) = 
plateau expected in the large length L-limit, actually 
emerges at (M) = 2/L. Fig. 3 illustrates this observation 
for L = 24, 20 and 16 and should emphasize the essential 
role of PBC in all our subsequent numerical checks. 

By construction, the bosonization approach can be 
straightforwardly extended to the case in which more 
than one frequency is present in e n . It turns out that 
whenever condition (|^) is satisfied for each frequency, a 
magnetization plateau shows up. This extension is simple 
since each operator (corresponding to each frequency) is 
commensurate separately, and hence the different pertur- 
bations can be treated separately. Of course the situation 
changes in the case of a dense multi-frequency spectrum 
(such as in the Fibonacci potential that we discuss be- 
low) and a more complete analysis has to be carried out 

& 

In order to test the reliability of these predictions, we 
have analyzed numerically both XX and XXZ chains 
using double-frequency couplings. In Fig. 4 we display 
the magnetization curves obtained for ui\ = 5/8, u>2 = 
7/8 with amplitudes Si = 0.2 and 6% = 0.3 respectively. 
As we mentioned before, the rather robust plateaux 
emerging at (M) =1/4 and 3/4 not only confirm the cor- 
rectness of our extended bosonization prediction, but also 
pave the way to an alternative description of the mas- 
sive spin excitations observed in NH4CUCI3; a quasi-one- 



dimensional S 



compound which is attracting both 



theoretical and experimental attention and whose magne- 
tization behavior yet remains unexplained |^,^,^5| . Of 
course, a more realistic description of this material must 
start from the microscopic structure observed by X-ray 
spectroscopy, which points to a two- leg zigzag ladder. 
Still, it is quite encouraging to obtain a magnetization 
curve qualitatively similar with the simple two-frequency 
model we considered here. 

In studying irrational frequencies or other quasi- 
periodic modulations, it is natural to check our analysis 
with the prototype Fibonacci sequence, a coupling array 
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J A = + 8), Jb = J(l ^ <S) generated by iterating the 
substitution rules B —> A and A -> AS [f|f||ll]] . Here we 
discuss the general XX Z situation (a related model has 
been studied in |B6[ ) , and compare our results with the 
already well known magnetization curve of the XX case 
[§,[5). Before continuing with the bosonization approach, 
we pause to discuss the strong coupling regimes of this 
system (8 — > ±1) in the context of a simple decimation 
procedure (see e.g. |^7j]). To evaluate the magnetization 
of the widest plateaux, there are two different cases to 
consider, according to 5 ~ —1, i.e. Jb 3> Jai an d the 
opposite situation for 8 ~ 1. 

Starting from saturation, in the first case the mag- 
netic field is lowered until it reaches the value h c ~ Jb 
at which the type-i? bonds experience a transition from 
the state of maximum polarization to the singlet state. 
The magnetization at this plateau is then obtained by 
decimating the B bonds. This simply yields (M) — 
1 — 2Nb/{Nb + Na), where Na,b denotes the number 
of bonds of type A and B respectively. For a large iter- 
ation number of the rules referred to above (L — > oo ) , 
Na/Nb approaches the golden mean 7 = (1 + n/5)/2 and 
therefore we find M l = (7 - l)/(7 + 1) ~ 0.236068. 

In the second case, J a Jb, we have to distinguish 
two different unit cells since type- A bonds can appear 
either in pairs (forming trimers) or isolated (forming 
dimers). It can be readily checked that when lowering 
the magnetic field from saturation the first spins to be 
decimated correspond to those forming trimers. We then 
find two plateaux at M 2 = 1 - 2/7 3 ~ 0.527864 (after 
decimating trimers) and, alike the case S ~ 1 , at M\ 
(after decimating the remaining dimers). Since the dec- 
imation procedure applies for generic XX Z chains p7[ , 
we conclude that the emergence of these strong coupling 
plateaux is a generic feature, at least with an antiferro- 
magnetic anisotropy parameter < A < 1. 

For intermediate regimes < \8\ < 1 the magnetiza- 
tion curve has a much richer structure which can be easily 
understood from our bosonization analysis in a multi- 
frequency case. Evidently, the self similar hierarchy of 
frequencies resulting from the Fourier transform of the 
Fibonacci exchanges (see e.g. along with the quan- 
tization condition studied so far, enables to reconstruct 
the whole spin gap structure of the Fibonacci chain, at 
least for | S | <C 1 and as long as the operator responsible 
for the plateaux is relevant. Actually, it turns out that 
one can obtain a fair approximation to the latter showing 
the most important plateaux by keeping just a few num- 
ber of main spectrum frequencies, even beyond weak cou- 
pling regimes. Interestingly, when the constraint (|^) is 
applied to the dominant Fibonacci frequencies 0^1 = 1/7 
and loi = 2(1 — 1/7), it yields precisely the M\ and 
M% plateaux arising from the strong coupling decimation. 
These results are shown in Fig. 5 where both single and 
double-frequency approximations are displayed. More- 
over, using the Fourier spectrum to set 81/ 8\ ~ 0.36237, 
the Fibonacci gap widths exhibit a remarkable agreement 



with the scaling exponents 1/(2 — d) referred to above. 
This can be observed in the inset over more than three 
decades in 8. In turn, this supports the observation that 
all gap widths of the XX Fibonacci chain scale simulta- 
neously with 8, as for A = the compactification radii 
comprehended in the scaling exponents are independent 
of the magnetization [R = l/(2y / 7r)] (see e.g. p"8||). 

Finally, we point out that for the Fibonacci chain 
(which has a dense Fourier spectrum), it was shown 
by solving the one-loop RG equations at zero magnetic 
field (jll| (see also |Q) that the critical value of R at 
which the Kostcrlitz-Thoulcss transition occurs is moved 
to R c = 1/V47T, in contrast to R c = 1/V87T for the sin- 
gle, or more generally non-dense frequency case. Hence, 
for < A < 1 and arbitrary magnetization, the op- 
erator responsible for spin gap openings at each of the 
Fibonacci frequencies is relevant. Therefore, the magne- 
tization curve of the XXZ case with antiferromagnetic 
anisotropy is of the same form as that of the XX situa- 
tion (see Fig. 5), though with different plateaux widths. 

To summarize, we have studied the interplay between 
quasi-periodic exchanges and uniform magnetic fields in 
strongly correlated antiferromagnetic chains using both 
bosonization and numerical techniques. The former were 
tested and complemented by the latter in a variety of non- 
perturbative scenarios. Our calculations suggest the pos- 
sibility to observe rather stable magnetization plateaux 
[Eq. (g) ] on artificially grown arrays of quantum dots ^ 
according to the (controlled) spatial distribution of their 
exchange integrals. We trust this work will convey an in- 
teresting motivation for further experimental studies in 
these material technologies. 
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FIG. 1. Critical fields (in bold lines) of the XX chain with 
uj = 13/21 and 10 spins. At 5 — 0.9 these fields conform the 
standard magnetization curve displayed in the inset. Vertical 
lines denote regions of massless spin excitations where the 
magnetization increases continuously with h. Empty zones 
in ascending order represent plateaux appearing at {M) = 
1/21, 5/21 [central region expected by Eq. (6)], 9/21, 11/21 
and 15/21 before saturation (top zone). 




FIG. 2. Magnetization contours of the isotropic (A = 1) 
XX Z chain for ui — 5/8 . Solid, dashed and dotted lines 
stand respectively for the critical fields of L = 24, 16 and 
8. The middle empty region corresponds to the (M) = 1/4 
plateau expected by Eq. (6). The inset shows one of the typ- 
ical magnetization curves upon setting 5 = 0.5 . 
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FIG. 3. Magnetization curves of a dimerized Heisenberg 
chain (a> = 1/2 with 5 = 0.4), using open boundary conditions 
displaying a plateau at (M) = 2/L . Solid, dashed and dotted 
lines denote respectively the results of L = 24, 20 and 16. 



FIG. 5. Fibonacci magnetization for 8 = 0.5 and 
F(25) = 75025 XX spins (bold line). Here, F(25) denotes the 
25 th Fibonacci number defined as F(n + 1) = F(n) + F(n-1) 
with F(l) = F(2) = 1. Dotted and solid lines stand respec- 
tively for the single and double frequency approximants to the 
Mi and M2 plateaux referred to in the text. The slopes com- 
pared in the inset show the gap width of M2 and Mi, in as- 
cending order. Bold lines denote the Fibonacci gaps whereas 
solid lines represent the double approximant widths. 




FIG. 4. Double frequency magnetization curve of the 
isotropic (A = 1) XXZ chain for u>\ = 5/8, UJ2 = 7/8 with 
amplitudes 81 = 0.2 and 82 = 0.3 . Solid and dashed lines de- 
note respectively the magnetizations of L — 24 and 16. The 
inset shows the magnetization curve of the corresponding XX 
chain with 10 5 spins. 
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